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Abstract 

Standard electroweak baryogenesis in the context of a first order phase transition 
is effective in generating the baryon asymmetry of the universe if the broken phase 
bubbles expand at subsonic speed, so that CP asymmetric currents can diffuse in front 
of the wall. Here we present a new mechanism for electroweak baryogenesis which 
operates for supersonic bubble walls. It relies on the formation of small bubbles of the 
symmetric phase behind the bubble wall, in the broken phase, due to the heating of 
the plasma as the wall passes by. We apply the mechanism to a model in which the 
Higgs field is coupled to several singlets, and find that enough baryon asymmetry is 
generated for reasonable values of the parameter space. 
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1 Introduction 

The most interesting aspects of a first order electroweak (EW) phase transition are that 
it may explain the observed baryon asymmetry of the universe through EW baryogenesis 
[TH3], and that it generates a stochastic background of gravitational waves (GWs) possibly 
detectable by next generation space-based interferometers [1H9]. The dynamics of the growth 
of the broken phase bubbles, and in particular the bubble wall velocity, play a fundamental 
role in both aspects. Fast walls are needed for sizable GW signals (the amplitude of the 
signal depending, among others, on the cube of the wall velocity v^, (TDHIS]), while slow 
walls are needed for EW baryogenesis: only subsonic walls (or a compression wave) lead to 
efficient diffusion of particle asymmetries in front of the wall |3j. 

An explicit computation of the bubble wall velocity has been done in the Standard Model 
(SM) [I3l and in the Minimal Supersymmetric Standard Model (MSSM) [14j, with the result 
that in the MSSM the bubbles are quite slow (f^ ~ 0.05 — 0.1), whereas in the SM the wall 
velocity is significantly larger, but still subsonic {vw ~ 0.35 — 0.45). However, for the SM 
the EW phase transition is not first order but rather a smooth cross-over [15] for values of 
the Higgs mass above the LEP bound rrih > 114.4 GeV [16], and MSSM EW baryogenesis is 
now strongly disfavoured by LEP and Tevatron data [17] and may be even more challenged 
by the LHC in the near future. Yet, there are other extensions of the SM where the three 
Sakharov conditions for baryogenesis are naturally satisfied during the EW phase transition, 
such as for example composite Higgs models [18], or extensions of the MSSM with a stronger 
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EW phase transition [1^1120] • Qualitative arguments show that in many of these extensions 
that lead to a first order EW phase transition, the wall velocity could be rather large, since 
it increases with the strength of the phase transition [21]. This happens typically in the 
case of non-supersymmetric extensions, where few degrees of freedom are added to the SM 
below or close to the TeV scale, but also in the case of supersymmetric extensions, if the first 
order phase transition is strong enough. Therefore, in these scenarios it might happen that 
the three Sakharov conditions are satisfied, yet standard EW baryogenesis cannot operate 
because the bubble wall velocity is large enough to suppress the diffusion of the particle 
asymmetries ahead of the bubble wall. 

In this perspective, we present here a new EW baryogenesis mechanism which operates 
in the case of supersonic bubble walls. We focus in particular on bubble expansion as a 
detonation. The mechanism is based on the fact that the expansion of the bubble as a 
detonation causes heating of the broken phase plasma in a small region behind the bubble 
wall [22j. If the temperature in this region exceeds the critical temperature Tc, bubbles 
of the symmetric phase can nucleate and grow inside the broken phase bubble. In these 
"symmetric" bubbles sphalerons are reactivated, due to the restoration of the EW symmetry, 
and baryon number is not conserved. The expanding symmetric bubbles move along the 
detonation wave towards the centre of the broken phase bubble, and at some point they 
enter the region where the heating is not very effective and the temperature is lower than 
Tc. In this region, the expansion of the symmetric bubbles is not any longer energetically 
favourable, and the bubbles start to shrink: the plasma then flows from the symmetric to 
the broken phase across the symmetric bubbles wall, and the baryon asymmetry generated 
inside the symmetric bubbles is transferred to the broken phase. We show that this process 
happens at a sufficiently low velocity to guarantee the diffusion of the asymmetry. For 
this mechanism to work, one further has to make sure that a large enough volume of the 
plasma goes through the symmetric phase, to explain the observed baryon asymmetry of the 
universe. We apply the mechanism to a model in which the Higgs field is coupled to several 
singlets, and analyse for which values of the parameters in the model (the universal coupling, 
the wall velocity, the Higgs mass) the filling factor is large enough and EW baryogenesis can 
occur for supersonic broken phase bubble walls. 

Having significant baryogenesis in models where the wall velocity is supersonic could 
accommodate a sizable GW signal with effective baryogenesis. Clearly, scenarios which can 
explain the baryon asymmetry of the universe and at the same time be tested with GW 
detectors are potentially very interesting (see for example [23]). Note that this can happen 
also in the case of subsonic wall velocities: in [2l] it has been pointed out that the relevant 
velocity for transport, being the velocity of the plasma as the wall sweeps through it, may 
be in general significantly lower than the wall velocity (although a more detailed analysis is 
needed) . 

To describe the evolution of the broken phase bubbles, we use the hydrodynamic de- 
scription of the bubble expansion (based on [25] and first applied to primordial cosmological 
phase transitions in [25]), where the plasma around the bubble is in thermal equilibrium 
and the bubble wall reaches a steady state with constant velocity, while the center of the 
bubble is at rest. However, this picture cannot be applied to the symmetric bubbles, since 
they expand in an inhomogeneous and anisotropic background (the detonation wave of the 
supersonic broken phase bubbles). We therefore develop a qualitative description of the sym- 
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metric bubble expansion, which is rehable and captures the relevant aspects of the problem: 
a more detailed treatment is beyond the scope of this paper. 

The paper is organised as follows: in sections |2] and [3] we review the hydrodynamics of 
bubble growth, the behaviour of the plasma as the detonation wave passes through, and 
the heating of the plasma in the rarefaction wave behind the bubble wall. In section H] we 
demonstrate that nucleation of symmetric bubbles is possible, and in section |5] we analyse 
how they grow after nucleation. Then, we turn to baryogenesis: in section [6] we discuss in 
detail the way in which baryogenesis takes place and evaluate the filling factor, and in section 
[7] we analize a specific scenario as an example in which supersonic baryogenesis is actually 
possible. Finally, we conclude in section |8l Some of the relevant but involved calculations 
are left for the appendixes. 



2 Hydrodynamics of Bubble Growth 

In this section we introduce the hydrodynamic analysis of the combined wall-plasma system 
(for a recent review, see that will later be used in sections [3] and O The two 

basic assumptions leading to the hydrodynamic approximation are local thermal equilibrium 
in the plasma and energy-momentum conservation in the systemH. The energy-momentum 
tensor of the Higgs field is given by 
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with Vo(0) being the renormalized vacuum potential. The energy-momentum tensor of a 
plasma locally in thermal equilibrium is 

T^™ = wu,u,-g,,p, (2) 

where w and p are the plasma enthalpy and pressure, and the quantity = (7, 7v) is the 
four-velocity field of the plasma. A constant background from ([T]) contributes to the total 
pressure, and will be included in p from now on. The enthalpy w, the entropy density a and 
the energy density e are defined by (T being the temperature of the plasma) 

_rj.dp _ dp _rj.dp , 

In the cases where the bubble expands at a constant speed, energy-momentum conserva- 
tion d^Tfj^u = d^T'l^^ + d^TJ^]^^^ = reads in the wall frame (with the wall and fluid velocities 
aligned in the z direction) dzT^^ = dzT^^ = 0. Integrating these equations across the bubble 
wall and denoting the phases by subscripts + (symmetric phase) and — (broken phase) one 
obtains the matching equations [25|[26] 



^The system is governed by a set of coupled equations of motion for the plasma and the Higgs field. 
In particular, the plasma is very well described by a fluid with phase space density f{p,x,t) obeying a 
Boltzmann equation. Then, when the interaction rates in the plasma are fast, the collision integral in the 
Boltzmann equation forces the phase space density towards a form that minimizes it, namely local thermal 
equilibrium, and the equations of motion for the plasma can be replaced by these two assumptions. 
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v+v- = , — = — . (4) 

e+ — e_ V- e+ + p- 

In the general case, the pressures and energy densities in the symmetric and broken 
phases are obtained from the free energy J-" = —p (see ([3])). In a concrete model (with 
known J-") the temperature T/v at which the phase transition happens can be determined 
using the standard techniques |271|28] (see section H]), and the thermodynamic potentials can 
be calculated in the two phases. Still, there are three unknown quantities (f+, w_ and either 
T_|_ or T_) and only two equations so that up to this point all hydrodynamically viable 
solutions are parametrized by one parameter, usually chosen to be the wall velocity f^. 

In order to obtain the behaviour of the system far from the wall, energy-momentum 
conservation in the plasma Qi^T^J^™^ = is imposeclfl, giving rise to a set of hydrodynamic 
equations [25l|26l[22]. Since there is no intrinsic macroscopic length scale present in the 
system, the solutions to the hydrodynamic equations are self-similar, depending only on the 
combination ^ = r/t (where r denotes the radial coordinate of the bubble from its centre 
and t the time since nucleation), and the hydrodynamic equations read 
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1 - v\i) [ cl ^\ i ' 1 - v\^) w{i) ^ ^ 

where Cs is the speed of sound of the plasma and /i is the Lorentz transformed plasma fluid 
velocity 

^^{^,v) = ^^. (6) 

1 — 4 f 

The solutions to correspond to the velocity profile f(^) and enthalpy profile w{^) 
(or temperature profile T{C,)) of the plasma away from the bubble wall. Quite generally, 
there are three different types of solutions (see [22] for details): detonations, deflagrations 
and hybrid solutions. In detonations j26], the wall expands at supersonic velocities > Cg 
and the vacuum energy of the Higgs leads to a rarefaction wave behind the bubble, while 
the plasma in front of the wall is at rest (see Figure [T]). In this case, the wall velocity is 

= f + > f _ > Cs [2S], and the temperature outside the bubble wall corresponds to the 
temperature of the universe at the time of the phase transition T+ = T^r. In deflagrations [30] , 
the plasma is mostly affected by reflection of particles at the bubble wall and a compression 
wave builds up in front of the wall, leading to T+ > T/v. "Pure" deflagrations are subsonic, 
with the plasma behind the wall at rest and = > v^^ while the hybrid case occurs 
for supersonic deflagrations [21] where both effects (compression and rarefaction wave) are 
present. 

We are interested in the behaviour of a plasma volume element in the case of a detonation, 
since it will be used in the following of the analysis. In this case, a plasma volume element 
initially at rest at a distance tq from the center of the bubble will feel the arrival of the 
detonation wave at time t = ro/f^ and will be dragged by it, moving along the wave 
according to 



•^Away from the wall, T/f^^ just gives a constant background that docs not play any role in energy- 



momentum conservation. 
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Figure 1: The solutions for the plasma velocity v{^) (LEFT) and temperature T{^)/Tn (RIGHT) 
in the case of a detonation (for Vyj = 0.76 and a^v = 0.078). The dashed-red horizontal line 
corresponds to Tjv and the solid-blue horizontal line corresponds to Tc- 



and eventually coming back to rest (for ^ = Cs) at a distance r > ro from the center of 
the bubble, as shown in Figure [2] (the motion of plasma volume elements in the case of 
deflagrations and hybrids can be described in a similar way). 

It is possible to extract the ^ and ro dependence of a trajectory r obtained from ([7]), 
using the self-similarity of the detonation wave. From Figure O we obtain the relation 

Arc, = r(ro + Aro, ^ = Cs) - r{ro, ^ = Cs) = C{cs) Avq , (8) 

where the constant C{cs) > 1 does not depend on tq. Due to the self-similarity, ([8]) can be 
generalized to arbitrary values of ^ inside the detonation wave 

Ar^ = r(ro + Aro, - rir^, = C(0 Aro , (9) 

where the function C(^) is again independent of ro, and one has C(^) > 1 all along the 
detonation wave. It then follows that 

r(ro, = C(0 ro + r(0, ^ C(0 ro , (10) 

where r(0,^) can be safely neglectecQ with respect to ro. The result ffTU]) shows the explicit 
dependence of a trajectory r(ro, C,) with the initial distance ro, and will be used later in 
section |5l 

''if the distance of a plasma volume element to the center of the bubble is negligible compared to the 
final radius Rb oi the bubble (ro — > for that volume element), then its distance to the bubble's center as 
it moves along the detonation wave is r(0,^) ^ tq, generically for values of tq ^ O(10~'^Rb) — 0{Rb)- 
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Figure 2: Motion of plasma volume elements (black lines) as the detonation wave (in grey) sweeps 
through them, for different values of the initial distance rg to the center of the bubble. 



3 Plasma Heating from Bubble Expansion 

From Figured] it appears that the plasma temperature behind the wall (in the broken phase) 
is higher than the nucleation temperature T^r, as a consequence of the expansion of the 
bubble. In this section we will perform a general analysis of the plasma heating behind the 
bubble wall, focusing on the case of bubbles expanding as detonations. Instead of performing 
the analysis in the general case p = —J-", we consider for simplicity the Bag EoS, holding if 
the plasma is a relativistic perfect fluid: 

P+ = ^«+^+ e+ = «+^+ + ^ ' (11) 

P- = e_ = a_Ti , (12) 

where e = Vo(0) — Vo(0o) denotes the false- vacuum energy resulting from the Higgs potential, 
and a± relate directly to the number of relativistic degrees of freedom g* in the symmetric and 
broken phases {g^ = 30 a-t/vr^). Using (fTTj) and (fT2l) we can write the matching equations 
(H as 

_ 1 - (1 - 3a+)6 f + _ 3 + (1 - 3a+)6 

~ 3-3(l + a+)6 ~ ~ l + 3(l + a+)6 ' ^ ' 



where we defined 



a+^^ b^^. (14) 
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The quantity a, being the rate of vacuum to thermal energy, typically characterizes the 
strength of the phase transition (the larger a the stronger the phase transition). Finally, the 
two equations in (|T3|l can be combined to give 



(15) 



1 + 




so that there are two branches of solutions, corresponding to the ± signs in ( fTSl) . The 
presence of a compression wave in front of the wall (deflagrations and hybrid solutions) 
corresponds to the — branch in (ITSj) . In this case, the temperature just in front of the wall 
T+ is higher than the temperature of the universe T/v- Also, depending on the strength of 
the transition, it is possible for the temperature just behind the wall T_ to be larger than 
Tat. However, there is a hydrodynamic obstruction to having T_ > Tc [32]. This is clear 
since the average temperature in the wall is T < (otherwise the net pressure on the wall 
that drives the expansion vanishes), and for a compression wave one 

For detonations the situation is quite different. In this case there is no compression wave 
in front of the wall, so T+ = T/v and Vj^ = > V- >_ Cg (with Cg = 1/ \/3 for a relativistic 
perfect fluid). The quantity b = (a^Tj^) / (a^T^^) , that measures the heating of the plasma 
behind the wall, is obtained from f[T^ 

b= — , — (16) 

If Vw denotes the lowest possible detonation wall velocity, obtained for t>_ = Cs, since VwCs < 
VyjV- < 1 it follows that h is bounded by 

'^CgVw — 1 < 6 < ^Yj\^ 

?)CsVw{l + a^v) — (1 — 3q;7v) 1 + ^Q^n 

From the upper bound on h it is clear that the plasma behind the wall is heated up with 
respect to the plasma in front, T_ > Tjy (since a+ > a J). We then want to know if for 
detonations it is possible to achieve sufficient heating for r_ to become larger than Tc. For 
the Bag EoS one has 

3 V a+y 1-^ 



Then, using (11611 the condition T_ > Tc reads 



J_ / a+ _ \ ?,v^v^{l + un) - (1 - 3ajv) ^ ^ 

Oat \a_ / "iv^Vyj — 1 



(19) 



The region of the (oat, f^)-plane (for a fixed ratio a_/a+ = 0.85) for which T_ > Tc is 
verified is shown in Figure [31 Apart from the trivial bound ajq > ac (T]y < Tc is needed 
for the bubbles to have been nucleated) there is a lower bound on (dashed-black line in 



^If T- > Tc would have been possible, this would have been potentially dangerous for EW baryogenesis 
scenarios due to the reactivation of sphaleron processes in the broken phase behind the bubble wall. 
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Figure 3: Region in the {aN,Vw)-plane for which T_ > Tc (fT9]) . for a-/a+ = 0.85. The dashed-red 
hne corresponds to Vw = Vw (lower bound on Vw), and the lower bound on a^v (dashed-black line) 
corresponds to the occurence of the hydrodynamic obstruction in the deflagration case |32j . The 
different values of fh (RIGHT) correspond to different sizes of the region where T > Tc with respect 
to the size of the bubble (i20l) . 



Figure [3]) due to the hydrodynamic obstruction that occurs in the deflagration cas (S. We 
now consider the temperature profile of the plasma T(^) in the rarefaction wave, evolving 
from T_ = T(f^) down to T{cs) < T^. In the cases when T_ > the condition T(^c) = Tc 
will be satisfied by some in the interval Cg < ^c < v^,, and the size of the region T(^) > Tc 
with respect to the size of the bubble will simply be 

h = . (20) 

From Figure [3] it appears that the region where the heating can be sufficient to drive T_ 
above Tc becomes smaller as the strength of the phase transition grows {i.e. for increasing 
a^). This counterintuitive effect is due to the interplay between and a^/a+. As shown 
in Figure m for fixed a^v the region in the (a_/a+, f^t,)-plane where T_ overcomes Tc becomes 
larger for decreasing a_/a_|_. This behaviour is general and occurs also for stronger phase 
transitions. In a reafistic model, and a_/a+ will not be independent quantities, and 
these two opposite effects will combine, either to reduce the heating effect for stronger phase 
transitions, or to enhance it depending on the details of the model (see section [7]). 

Having established that for detonations it is possible to heat the plasma behind the wall 
above Tc, we now proceed to study the consequences of this heating effect: in particular, we 
investigate whether the system can go back to the symmetric phase locally in the regions 
where T(^) > %. 

^For aN below tliis bound, the bubble will always expand as a subsonic deflagration, no matter how small 
the friction on the wall (due to the plasma) is. Therefore, a detonation is never realized in this case [32] ■ 
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Figure 4: Region in the (a_/a+, fu;)-plane for which T_ > Tc (fT9]) . for = 0.1. As in Figure [3l 
the dashed-red line corresponds to Vw = (lower bound on v^)- 

4 Nucleation of Symmetric Bubbles. 
4.1 Bubble Nucleation in the Standard Case. 

Consider the behaviour with temperature of a potential V{(f),T) giving rise to a first order 
phase transition in the early universe (being the order parameter of the transition), such 
as the one depicted in Figure |5l For very high temperatures the only minimum is at = 
(Figure |5] (a)). As the universe expands and the temperature lowers, a new minimum 0^ 
develops away from the origin and eventually the system reaches a critical temperature Tc 
for which V^(0) = V^cps) (Figure [5] (b)). Then, for T < Tc the potential has a metastable 
minimum at = and a stable one at 7^ 0, separated by a potential barrier (Figure |5](c)) 
and it becomes possible to tunnel from = to the stable minimum at 0b. The nucleation 
probability per unit time and volume from the metastable minimum to the stable one is 
given at finite temperature by [28] 

r,r).r(|g^)"\xp(-^). (21, 

with 5*3 (T) being the tunneling three dimensional Euclidean action |27ll28]. Then, in the 
radiation dominated era the mean number of true (broken) vacuum bubbles iVsubbies nucle- 
ated in a volume V during a time interval t — tc (being tc the time at which T = Tc) would 
be simply 

iVBubbies= / T{t)V{t)dt= / I ^ V{T)T{T)dT, (22) 

Jtc Jt -L 

where the Hubble rate H{T) and the time-temperature relation during the radiation domi- 
nated era read 
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Figure 5: LEFT: Evolution of tlie potential V{(I),T) with temperature in the case of a first order 
phase transition in the early universe (see section |4?T]) . RIGHT: Local change in the potential for 
the case r_ > Tc (see section 14. 2p . 

Even though nucleation of broken phase bubbles becomes possible as soon as the tem- 
perature drops below the critical temperature T^, the actual nucleation temperature T^r at 
which the phase transition starts is defined as the one for which the mean number of bubbles 
nucleated inside a Hubble volume H~^{T) is of C(l) [33] 

^B..-(r„) = £ ^ (1^) exp (- dT . 1 . (24) 

and this in turn leads to the condition 

^ 142 . (25) 
4.2 Symmetric Bubble Nucleation. 

Consider now one of these nucleated bubbles of broken phase expanding in the plasma (from 
now on we will refer to these bubbles as "broken" bubbles). Recalling section 3, if the 
bubble expands as a detonation the plasma behind the bubble wall gets heated up, and if 
T_ > Tc (being T_ = T(w^) the temperature just behind the bubble wall) then for some 
region < ^ < i;^ of the detonation wave one has T{^) > Tc (with T(^c) = Tc). In 
this region the system is then in situation (d) of Figure [5] (RIGHT), and can tunnel back 
to the symmetric minimum, with the subsequent formation of bubbles of the symmetric 
phase inside the broken bubble (from now on we will refer to these bubbles as "symmetric" 
bubbles). Differently from the case discussed in the previous section, the available volume 
for the inverse tunneling and subsequent formation of the symmetric bubbles is not the entire 
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Hubble volume but just the region inside the broken bubble where T(^) > Tc, whose volume 
is Vf^. Recalling fl20|) and the fact that at a given time t after the nucleation of the broken 
bubble its radius will be Rb = v^t, (and we also define R^^ = C,ct), Vf^ can be written as 

Vf, = ^Rl - ^Rl = K - e.) ^A.{v^tfh{l-h). (26) 

The nucleation probability ( 12T|) is not constant over Vj^ (due to the temperature variation 
along the detonation wave), and so the mean number of nucleated symmetric bubbles inside 
a broken bubble is not defined by ( !22l) . Instead, the number of symmetric bubbles nucleated 
inside the broken bubble per unit time is obtained by integrating ( 12T]) over the region of 
volume V/j, 



l-tiB 

Ni{t)=ATT T{T{R))R^dR, (27) 

which then gives 

Ni{t)=47ct' r(T(0)e'^^e = 4 7rtM r (T) ^ e(T)' . (28) 

The total number of symmetric bubbles nucleated inside the broken bubble up to a time 
r (since the nucleation of the broken bubble^]) A/'Bybbies(''') J^^t found by integrating fj28|) 
over time 

A/'lubbiesM= f dtNi{t). (29) 
Recalling ( 12^ and ( 125|) . (129|) can be approximately rewritten as 

rs _A.o Sl{T^)\ f^' ^ (T) 



A/Bubbies(r) ^ jp-, exp (^142 - J ^ e(T)^ ciT , (30) 

where S^{T) is the tunneling action from the broken to the symmetric phase. Since r < (3^^, 
where is the duration of the phase transition {P/H ~ 10 — 1000 [35]), then from (l30l) the 
value of 5*^ (T_)/T_ needed to nucleate at least one symmetric bubble per broken bubble is 
bounded by 

-^3^^-) < 142 -4 log f 4^125, (31) 



where we have neglected the contribution from the integral in ( 130|) : this would in principle 
further suppress the above value for S'f (T_)/T_ due to the fact that the integrand is strongly 
peaked at T — )■ T_, but this would enter in ( 13T|) only logarithmically. The above equation 
shows that, due to the volume suppression, it is generically more difficult to tunnel to the 
symmetric phase. 



^Actually, since symmetric bubbles cannot nucleate until the broken bubble reaches a stationary state 
and the detonation wave is established, r is really the time since the broken bubble has reached a stationary 
state, and not the time since nucleation. However, both are the same for all practical purposes, since the 
stationary state is reached extremely fast compared to the duration of the phase transition. 
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A symmetric bubble can in principle be nucleated anywhere in the volume Vf^. However, 
the situation is different from the case of broken bubble nucleation, which is a true random 
process since the nucleation probability is the same everywhere. Here, the nucleation proba- 
bility is much higher close to ^ = f ^ than close to ^ = ^c, because one has T{vw) > T{C,c) and 
the nucleation probability r(T) varies exponentially with T ( 12T1) . The average nucleation 
value for ^ is then very close to Vw 

Most of the symmetric bubbles are then nucleated very close to the broken bubble wall. 



5 Growth of Symmetric Bubbles. 

Let us assume that a symmetric bubble nucleates at S,s — (at a given time t^^^^ after 
the nucleation of the broken bubble) and subsequently expands. Due to the inhomogeneity 
and anisotropy of the background in which the symmetric bubbles grow, the expansion will 
render them non-spherical. As a consequence, the symmetric bubble expansion cannot be 
described by a self-similar solution in a way analogous to the case of the broken bubble 
(recall section [2]). In fact, a rigorous solution of the system including the symmetric bubble 
expanding inside the broken bubble would require solving d^T^^ = over the whole system 
(with boundary conditions on the symmetric and broken bubble walls), a very difficult 
problem both analytically and numerically. We expect however that a good approximation 
is to consider the broken bubble as a background on which the symmetric bubble expands, 
neglecting the back-reaction on this background due to the symmetric bubble expansion. 
This is certainly the case in regions of the detonation wave of the broken bubble far from 
the symmetric bubble. Moreover, the symmetric bubble expansion turns out to be very slow 
(as we will see later in this section), and so it is reasonable to expect that the perturbation 
of the background is small. This is indeed similar to the case of broken bubbles expanding 
as very slow deflagrations, where the perturbation of the temperature background due to 
the bubble expansion is found to be very small [22] (not surprisingly, since the perturbation 
has to vanish in the limit — > 0). In the following we therefore consider that back- reaction 
is small and we neglect it, since this is justified even when a large amount of symmetric 
bubbles are nucleatecj^. 

Even though for T > the free energy of the system in the symmetric phase J^{(f) = 0, T) 
is smaller than the one in the broken phase J^{(j)B,T) (and inverse tunneling does indeed 
occur), the latent heat of the inverse transition is negative, due to the fact that for the 
symmetric bubble the broken and symmetric phases are interchanged (see Appendix A for 
details). Then, it is clear that the expansion of the symmetric bubbles cannot be fueled 
by the liberation of latent heat as the bubble expands (this is analogous to an endothermic 
reaction in the context of slow combustion [25|), and in this case one might be doubtful that 
the expansion of the symmetric bubbles actually takes place. 

^We also disregard here the perturbation on the background due to possible collisions between the sym- 
metric bubbles. These collisions are essentially absent when few symmetric bubbles are nucleated inside each 
broken bubble, but some collisions may occur if the number of symmetric bubbles is very large. 
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However, there are two facts that question this last conclusion. First, there is an increase 
in the thermal energy of the system as the symmetric bubbles expand (due to the increase 
in the number of relativistic degrees of freedom of the plasma when it gets converted to the 
symmetric phase) , which contributes to the pressure difference on the symmetric bubble wall 
driving the expansion. Second, since the heating of the plasma due to the expansion of the 
broken bubble keeps the temperature of the background in which the symmetric bubbles 
expand above for < ^ < Vyj, the temperature of the plasma in this region is effectively 
maintained despite the endothermic nature (negative latent heat) of the symmetric bubble 
expansion, making the expansion possible. In this scenario the expansion velocity of the 
symmetric bubble with respect to its background can be estimated qualitatively to be 

5v{i) = ^-l. (33) 

This can be understood as follows (as argued in [1]): as the symmetric bubble expands, 
the endothermic nature of the expansion will cause the temperature close to the symmetric 
bubble wall to drop momentarily. If the temperature would drop down to Tc, this would 
prevent further expansion [32]. However, the background temperature T{S,) close to the wall 
gets re-established by heat transport (since the expansion of the broken bubble keeps the 
background constant) and the symmetric bubble expansion continues. Since the rate of heat 
transport is roughly proportional to {T{^) — T(.)/Tc, the velocity of symmetric bubble growth 
will then be of the order of (T(^) — Tc)/Tc if the temperature difference is small^ |i|l25|l32]. 

Symmetric bubbles nucleate at rest with respect to the plasma in which nucleation occurs, 
and move along with it. The evolution of a symmetric bubble is shown in Figures E] and |Hi 
Recalling ([7]), the plasma element in which the symmetric bubble will be nucleated moves 
along the detonation wave of the broken bubble according to 

^=^i^ = ^c/t) , (34) 

where v{^) is the velocity profile of the detonation wave and rc{ro,C,) is the distance between 
the plasma element and the center of the broken bubblj^ (being tq its initial value, when 
the plasma element is at rest before the arrival of the detonation wave - see section |2]). 
The plasma volume element in which the symmetric bubble is nucleated passes through the 
detonation wave following a trajectory like those depicted in Figure |2] (solid-black lines), 
moving initially through the region where T{^) > Tc and entering the region with T(^) < Tc 
at a later time (see Figure [6]). 

Then, focusing on the expansion of the symmetric bubble along the broken bubble's 
radial direction (see Figure [7]) , the point of the symmetric bubble wall closest to the broken 
bubble center and the point of the symmetric bubble wall closest to the broken bubble 
wall Tout evolve as 

^Note that the situation described here is reversed from the one of A , which deals with the expansion 
of usual broken bubbles and analyses the case of ambient temperature T < Tc- In that case, the liberation 
of the (positive) latent heat causes the plasma right in front of the wall to be heated up to Tc, momentarily 
preventing further expansion. The liberated heat has then to be carried away into the symmetric phase for 
the expansion to continue. 

'^'^It is also effectively the distance between the center of the symmetric bubble and the center of the broken 
bubble (see Figure [7]) 
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Figure 6: Detonation wave of the broken bubble (in grey), and the symmetric bubble nucleated 
inside (in red). The solid-black line corresponds to the evolution of Tc (j34|) (for Vw = 0.76 and 
UN = 0.078). 



dt 1 _ ^; (lia) fa (Im) dt 1 + ^; (Imit) fa (I^) ' ^ ^ 

with 5v given by (l33l) . As seen from Figure [HI while the symmetric bubble is in the region 
^ > ^c, where T(^) > Tc, the plasma flows from the outside to the inside of the symmetric 
bubble (dashed-black line in Figure |8]), since Sv{^) > 0. However, once the symmetric 
bubble enters the region where T{^) < T^, the velocity 6v{^) becomes negative and the 
plasma inside the symmetric bubble, as seen from the sjTiimetric bubble wall, moves towards 
the wall instead of away from it. Then, the plasma flows from the inside to the outside 
of the symmetric bubble (dashed-black line in Figure [S]) and the symmetric bubble starts 
contracting, shrinking progressively and eventually disappearing. 

Given that a symmetric bubble expands to reach its maximal size for ^ = and then 
shrinks back to nothing, only plasma volume elements whose trajectories rc{ro,C,) are close 
enough to the one of the plasma element in which the symmetric bubble nucleates, rc(ro, 
will actually feel the presence of the symmetric bubble by going through it (dashed-black 
line in Figure [8]). The maximum and minimum value of fg are, respectively 

^o^ax = '^o + Sr^o^ ^o„,i„ = '^o - Sr^o^ , (36) 

where 6r^^ and 6r^^ are obtained from the conditions (see Figure E]) 

ront{ro, Q = rdro + Sr^Q\Q rm{ro, Q = r^ivo - 5rJ^\ ^c) • (37) 
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Figure 7: Sketch of the expansion of a symmetric bubble inside a broken bubble, and definition of 
rc, Tin and rout- 



Getting Sr^ and 5rQ as a function of vq and ^ is rather involved; we present the calculation 
in Appendix B. 

So far we have focused only on the evolution of the symmetric bubble along the broken 
bubble's radial direction, disregarding the expansion of the symmetric bubble along the 
angular directions. The full three-dimensional problem, together with the departure from 
the initial spherical shape of the symmetric bubble as it evolves, are difficult matters to study 
(see Appendix B). However, given spherical symmetry it is a reasonable approximation to 
consider that the volume of plasma that goes through a symmetric bubble is roughly given 
by dEl]) 

Vsiro) ^ (5r^) + 5rf )3 = {6rof . (38) 

This last result will be used in the next section, where we analyse the possibility of achieving 
enough baryogenesis through the nucleation of symmetric bubbles. 

6 Fast Bubble Expansion and EW Baryogenesis. 

The standard EW baryogenesis mechanism is based on the interaction between the wall of 
the expanding bubbles and the plasma in front of it [IHS], which leads to a CP-asymmetric 
reflection on the wall of certain particle species. These particle asymmetries are subsequently 
diffused into the plasma in the symmetric phase just in front of the bubble wall [3], where 
sphalerons are active and capable of converting the CP asymmetry into a net baryon number. 
The generated baryon number is then carried into the broken phase (where it stays frozen) 
as the wall passes by. 

The diffusion process plays a key role in the generation of the baryon asymmetry, since 
it connects the non-equilibrium CP violating physics (occurring within the bubble wall) to 
the sphaleron baryon number violating processes (occurring in the symmetric phase out- 
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Figure 8: Detail of Figure [6l The symmetric bubble is depicted in red, its upper boundary 
corresponding to the evolution of rout (|35p and its lower boundary corresponding to the evolution 
of Tin (j35p . Plasma volume elements with trajectories rc(fo,^) contained between the upper and 
lower solid-black lines go through the symmetric bubble (dashed-black line). 



side the bubble). Therefore, in order for the whole mechanism to be able to generate the 
baryon asymmetry of the universe, the diffusion process has to be effective, meaning that 
the timescale for the diffusion of the asymmetries into the symmetric phase has to be smaller 
than the time the wall takes to sweep through the plasma just in front of it (suppressing the 
sphaleron processes as it passes through). This puts an upper bound on the relative velocity 
between the bubble wall and the plasma in the symmetric phase v^. < D/Ly^ ~ 0.15 — 0.3 
(with D a diffusion constant and the wall thickness) [3]. Moreover, effective diffusion 
cannot take place if the relative velocity Vr is greater than the speed of sound of the plasma 
Cg. This seems to prevent the generation of the baryon asymmetry of the universe via EW 
baryogenesis (even if all three Sakharov conditions are present at the EW scale) for scenarios 
with relatively fast bubble walls, and in particular when the bubbles expand as detonations, 
case in which the bubble walls are supersonic!"]. 

Nevertheless, even if diffusion in front of the broken phase bubble wall is ineffective, and 
the EW baryogenesis mechanism previously described does not work, it may still be possible 
to generate the baryon asymmetry when the bubbles expand as detonations. Recalling the 
discussion of section [3|, the heating of the plasma behind the broken bubble wall caused by 

^^As explained in the introduction, while this is not the case for the MSSM, due to the phase transition 
being at most weakly first order and the large amount of friction on the wall |14| , it may indeed be the case 
for other extensions of the Standard Model with a stronger EW phase transition, such as the NMSSM [20] 
or the Standard Model extended with a singlet scalar [T8 l [33 ll38H42] . 
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the detonation wave may rise the temperature of a region behind the wall above and small 
bubbles of the symmetric phase will be nucleated in this region (see section in whose 
interior sphalerons processes are unsupressed (even though the symmetric bubble maximal 
size is always much smaller than the size of the broken bubble, as seen from Figure [6l it is still 
larger than the characteristic size of a sphaleron configuration by many orders of magnitude). 
Since this small symmetric bubbles are nucleated at rest in the plasma, they move along with 
the plasma volume element in which they are nucleated, and therefore they go through the 
detonation wave. While they are in the region where T > T^, the symmetric bubbles will 
slowly expand, with the plasma flowing across their wall from the broken to the symmetric 
phase. However, the symmetric bubbles reach a maximum size and then start shrinking 
as they enter the region where T < Tc, getting smaller until eventually disappearing (see 
section [5]). During the shrinking process, the plasma flows across the symmetric bubble wall 
from the symmetric to the broken phase. 



0,1 
0,05 

8v(^) 

-0,05 



-0,1 

0,55 C, 0,6 0,65 0.7 0,75 

Figure 9: 5v{S,) for a nucleated symmetric bubble (for = 0.76 and aj\f = 0.078). The red line 
corresponds to ^ < so that 5w(^) becomes negative, the plasma flows from the symmetric phase 
to the broken phase and baryogenesis can take place. We see that the range of velocities involved 
in this case is \Sv\ ~ 0.01 — 0.05 <C Cg. 

During the symmetric bubble expansion, EW baryogenesis is not possible, since the 
plasma flow across the bubble wall goes from the broken to the symmetric phase and there is 
no way any baryon asymmetry generated by sphalerons can successfully be transported into 
the broken phase before it is washed out. On the other hand, during the symmetric bubble 
contraction, the plasma flows across the wall from the symmetric to the broken phase (just 
as for the expansion of the broken bubble) and any potential baryon asymmetry generated 
by sphalerons in the symmetric phase is automatically transported into the broken phase as 
the symmetric bubble wall passes through. Moreover, as opposed to the broken bubble (for 
which diffusion is suppressed due to > Cg), for the symmetric bubbles the relative velocity 
between the wall and the plasma in the symmetric phase as they contract is \6v\ -C (see 
( l33|) and Figure |9]), and so the diffusion of the particle asymmetries into the symmetric phase 
is indeed effective. 
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Let us now suppose that the theory beyond the Standard Model responsible for the first 
order EW phase transition has all the necessary ingredients for successfully generating the 
observed baryon asymmetry of the universe via EW baryogenesis (sufficient CP violation 
and phase transition sufficiently strong) , but that the expansion of the broken phase bubbles 
is supersonic and proceeding via detonations. The EW baryogenesis mechanism can then 
still take place close to the surface of the symmetric bubbles. Since the symmetric bubble 
wall is expected to be in all respects similar to the broken bubble wall, the baryogenesis 
mechanism at the microscopic level is the same as in the standard case; therefore, it can 
operate due to the small expansion velocity of the symmetric bubbles. 

Still, the EW baryogenesis scenario presented here needs to address a potential "filling 
factor" suppression: this is due to the fact that, depending on the number of nucleated 
symmetric bubbles, it is not guaranteed that all the plasma inside the broken bubble will 
have undergone EW baryogenesis at one time or another by going through a symmetric 
bubble. If we define ub to be the local baryon number density generated during the process 
of EW baryogenesis and ub to be the average baryon number density in the universe after 
the EW phase transition has completed, then in the present baryogenesis scenario they are 
related by 

^Sym 

nB = riB , (39) 

where is the volume of a broken bubble and V^^^ is the total volume of plasma within the 
broken bubble having undergone EW baryogenesis by going through the symmetric bubbles. 
We then define the filling factor T = V^^'^/V'^. Since and y^y™ ultimately depends on 
the phase transition parameters (such as or T/v) and on the wall velocity w^, ( 139|1 can be 
written a^ 

^ = T(T„t;J. (40) 
Ub 

In order to get a positive filling factor, a necessary condition is that T_ > Tc (see Figure 
|3]). However, this is not a sufficient condition: one still needs to nucleate at least one 
symmetric bubble inside each broken bubble during the phase transition, assuming that the 
broken bubble was nucleated at the beginning of the phase transition (this assumption is 
well justified since those bubbles are in any case the ones that eventually cover most of the 
horizon volume at the end of the phase transition). This provides a lower bound on T_ (as 
a function of f^) for T > 0. From (1301) . this condition reads 

A/'ifubbics(r') = i- (41) 

Moreover, since the volume of plasma going through one symmetric bubble is very small 
compared to the size of the broken bubble at any time, Vsif) <C (see ( l38l) and the 
discussion in section |5]), it is only possible to obtain a sizable value of the filling factor T < 1 
if a large number of symmetric bubbles are nucleated inside each broken bubble. If only a 
few symmetric bubbles are nucleated per broken bubble, then T ^ 1 and hb will be highly 

""^^ While not shown exphcitly, T depends on other phase transition parameters apart from Tc- 
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suppresseco. Note that given the present uncertainties in the computation of in EW 
baryogenesis scenarios |36l[37], a reasonable requirement for successful EW baryogenesis is 
T > 0.1-0.5. 

Let us now turn to the actual computation of the filling factor T = V^^^/V^. Since the 
final radius of a typical broken bubble is Rb — (3~^, we have 



(42) 



On the other hand, the total volume of plasma within a broken bubble having undergone 
EW baryogenesis can be written as 



Vsit)Niit)dt^ 



(43) 



where Vs{t) = (5ro(t))'^ is the volume of plasma that goes through a symmetric bubble, 



obtained from (138|) and (IB.Sp . and N^(t) is the number of symmetric bubbles nucleated 
inside a broken one per unit time, given by (128|) . The lower integration limit tmm, defined 
by the condition 



1, 



(44) 



represents the fact that one needs to nucleate at least one symmetric bubble during the 
time interval Sro{t)/viu for Vs{t) to be nonzero (and EW baryogenesis to happen during that 
time interval). Moreover, the more symmetric bubbles are nucleated during a time interval 
6ro(t)/vuj, the less volume is left for further symmetric bubbles to nucleate, and so the value 
of the integrand in (H3|) cannot be arbitrarily high, but is bounded by 



Vs{t)Ni{t) 



Sro{t) 



rev, 



t - 



Sroit) 



(45) 



where the r.h.s of (H5l) corresponds to the total volume a broken bubble sweeps in a time 
interval (5ro(t)/f^, and no more symmetric bubbles may be nucleated in that volume. Since 
5ro(t) oc t and A^b(^) ^ if the bound ([15]) is reached for t = r < then (jl3]) should be 
replaced by 



ySym 



Vs{t)Ni{t)dt + -7ivl [/3- 



T 



(46) 



It is now clear that, for a given model that produces a first order EW phase transition, it 
is possible to obtain both V^^'^ and as a function of the parameters of the transition and 
the wall velocity v^,, using the results from sections [3l |4]and|5l In the next section we will 
study the efficiency of the proposed EW baryogenesis mechanism (including the evaluation 
of the filling factor T) in an explicit model. 



^■'This is similar to the case of baryogenesis from cosmic strings, for which the baryon asymmetry produced 
is volume-suppressed [34] with respect to ns due to the fact that the strings sweep just some small part of 
the volume of the universe as they move. 
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7 Supersonic Baryogenesis in a Specific Model. 



The analysis of the heating of the plasma due to the detonation wave (section [3]) can be 
performed without relying on any particular symmetry breaking scenario, simply using the 
hydrodynamic analysis from section [2] and specifying an equation of state for the plasma. 
On the other hand, the analysis of the nucleation and expansion of symmetric bubbles 
inside a broken bubble requires the knowledge of the free energy of the system J-", since the 
computation of the tunneling action can only be performed once the free energy has been 
specified. Consequently, in order to explore the viability of the EW baryogenesis mechanism 
here proposed, we concentrate on a specific scenario: the Standard Model extended by a set 
of singlet scalar fields S^. As explained below, for the sake of simplicity we take a particular 
limit of this model where the analysis can be performed analytically, and discuss afterwards 
the possible implications of a more general situation. 



7.1 Singlet Scalar Extension of the Standard Model. 



Extending the Standard Model by one or several singlet scalar fields coupled to the Higgs 
field yields a very attractive scenario in particular concerning EW baryogenesis, since this 
can easily give rise to a rather strong first order EW phase transition |33 [ [381440] . If one 
adds only one extra singlet field (either real or complex), the phase transition is sufficiently 
strong only if both the Higgs and the singlet change their VEVs during the phase transition 
(see for example [39lll0]). On the other hand, when the Higgs is the only field driving the 
phase transition (and for a Higgs mass above the LEP bound rrih > 114.4 GeV |16]) the 
presence of one extra singlet is not sufficient [41], and a relatively large number of singlets 
Ns is needed in order to achieve a strong phase transition 

The tunneling action 5*3 in a phase transition involving both the singlet field(s) and 
the Higgs field is difficult to obtain, and in general cannot be found analytically (see |13]). 
Therefore, for simplicity we consider here the particular case in which the Higgs is the only 
field changing VEV during the phase transition. This allows us to use a semi-analytical 
solution for S^. Then, in order to achieve a sufficiently strong first order phase transition, 
we have to add many singlet scalar fields (we take Ns = 8). We start from the following 
tree-level scalar potential 



(47) 



where for simplicity we assume that the mass term m?s-Sf is absent, and we take an universal 
coupling between the Higgs and the singlet scalars. 

Since in the above model ( H7|) the scalar singlet fields do not take a VEV, we can con- 
centrate on the potential along the Higgs field direction. In the background Higgs field 
configuration defined by {H^) = h/\/2, the effective potential at 1-loop in Landau gauge 
and MS renormalization scheme is 
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The subscript a = {Z,W,t, S^} denotes the gauge bosons and W"^), top quark 
and singlet scalar fields with Na = {3, 6, —12, Ns} (for the quantitative analysis, we take 
Ns = 8), while Ca = 5/6 (3/2) for gauge bosons (fermions and scalars). We are neglecting 
the contributions from all the quarks and leptons except the top quark, and also from the 
Higgs itself and the Goldstone bosons. The /i-dependent tree-level masses are 

Mlih) = Ch\ Ml{h) = \{g' + g'')h\ Ml{h) = \g'h\ M^{h) = \ylh\ (49) 

where g and g' are the SM gauge couplings and yt the top quark Yukawa coupling. The 
renormalization scale Q is chosen to be the top mass Q = Mt{v), and the parameters at that 
scale are fixed to recover (h) = v = 246 GeV. 

At a finite temperature T, the free energy of the Higgs field at 1-loop is given by 

T) = Vi_ioop(/.) + Y.^l^d^^" (l ± e-V^'^^^^WJ^) . (50) 

If we consider the high temperature limit T 3> Ma{h), we can perform an expansion in 
M^{h)/T^ in dSO]), obtaining 

aGbosons ctSfcrmions 

where g* = idb + (7/8)(7/) now includes also the contribution to the number of relativistic 
degrees of freedom of all the quark and leptons, as well as the Higgs and Goldstone bosons. 
Also, log(aa) = 5.41 (2.64) for bosons (fermions). 

Note that both in ( l50|) and ( ISTl) we are not including the contribution to the free energy 
of the Higgs field coming from the resummation of bosonic zero modes (the so-called "Ring" 
or "Daisy" contribution to J-"). This contribution is potentially important, specially at high 
temperature, and in a complete study of the phase transition it should be included; here we 
neglect it since we want to be able to perform an analytic study. Due precisely to the fact 
that we are ignoring this contribution, and that we are neglecting the contributions from the 
Higgs and Goldstone bosons in psj) and f lSOjl . we can write flSTl) as 

J^{h, T) = A{T) h'^ - B{T) + At(T) + (field independent terms) , (52) 

with 
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Arranging (IFIj) in the form (132]) makes the analysis much easier, as we now show. Eq. 
describes a system with a first order phase transition, the /i^-term being responsible for the 
potential barrier between minima. The extrema of fl52l) are: 




^ = h^ = VV^±\ -Trrh^- (53) 



The behaviour of the free energy (1321) with temperature is the following: at very high temper- 
atures the free energy has a parabolic-like shape, with h = being its only real extremum (a 
minimum), while h±. are complex. As the temperature lowers, it eventually reaches T = T^,, 
at which the term inside the square-root in (!53|) vanishes and the extrema h± become real 

Below T*, the free energy has two real minima [h = and hj^ ^ 0, with h = still being 
the absolute minimum) separated by a potential barrier (its maximum being at h^). As the 
temperature keeps lowering, it reaches T^, the temperature at which the two minima have 
equal free energy J^{h = 0,Tc) = J^{h^{Tc),Tc) 

52(T,) 32 4' ^ ^ 

For T < Tc the absolute minimum is h^, and it becomes possible to tunnel from h = to 
h^. In general, the computation of the tunneling action S3 has to be done numerically, but 
for the simple case of quartic potentials of a single scalar field 

V{(f)) = X(f)^ -b(f)^ + a(j)'^ A,a,6>0 (56) 

there exists a semi-analytic solution [H] for the three-dimensional Euclidean action S3, given 
by 

^^=A^ 81 (2-^)2 (/^i^ + /^2^' + /^3^') ^^-^ (57) 

where /3i = 8.2938, /^a = -5.5330, /^g = 0.8180. Since Eq. rewritten as in (132]), 

is precisely of the form (13^ (dropping field independent terms), we can use Eq. (137]) to 
compute the tunneling action. Then, the temperature T/v at which bubbles of the broken 
phase begin to nucleate is easily obtained from ( 124]) and ( 125]) . 

If the nucleated bubbles expand as detonations, the condition (fT9]) determines the region 
of parameter space where the temperature just behind the bubble wall T_ is larger than 
Tc. Note that now a_ and are obtained from J-' |22]- For the particular case we are 
analysing, once we fix the Higgs mass m/j, the strength of the first order phase transition 
is fully controlled by the value of the coupling ( between the Higgs and the scalar singlets 
(see Eq. (jS])). The larger is (, the stronger is the phase transition. We present our results 
letting ( and v^, as free parameters. In Figure [HJ] we show the region in the {(,Vw)-p\ane 
where T_ overcomes T^, for several values of rrih and bearing in mind that a phase transition 
sufficiently strong for viable EW baryogenesis needs Rc = h^{Tc)/Tc > 1. As previously 
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stressed, however, T_ > is a necessary condition for the supersonic EW baryogenesis 
mechanism to operate but not a sufficient one. 
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Figure 10: Allowed region for supersonic EW baryogenesis (in white) for the free energy (jSlh with 
nih = 125 GeV (up-left plot), rrih = 150 GeV (up-right plot), nih = 175 GeV (down-left plot) and 
nih = 200 GeV (down-right plot). In the yellow region, the stationary stable solutions are not 
detonations, but hybrids (see main text). 



For Tc < T_ < T* we can compute the tunnehng action from the broken to the symmetric 
minimum S^{T). This can be done by performing the shift h = — {h^ — h) = — 6h in 
Eq. (|52|) . which becomes then J^{h,T) = J^{h+,T) + J^{6h,T), with 
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T{5h,T) = XT{T){5h)^-[AXT{T)h+-B{T)]{5hf + 

+ [A{T) - 3 B{T) h+ + 6 XriT) hi] {Shf . (58) 

Again, this potential is of the form (!56|) . allowing us to use (!57|) to compute S'f (T). Once 
we know the Euclidean three-action, we can use the results of sections HI [5] and [6] to obtain 
the region in the (C,f^)-plane where symmetric bubbles are nucleated {J^iuhhiesif^'^) — 1) 
and thus the filling factor is non-zero T > (this region lies obviously inside the one for 
which T_ > Tc). We then compute the filling factor T in this region. The results are 
depicted in Figure [TOl which shows, for several values of the Higgs mass rrih, the region in 
the (C,f^)-plane where the proposed EW baryogenesis mechanism operates (in white). 

For values of below the dashed- red line in Figure [TU] (yellow region), no stable deto- 
nation solutions exist, since above the broken minimum ceases to exist. In this region, 
the solution most likely develops a compression wave in front of the wall and turns into a 
hybrid^ (or supersonic deflagration [31]), for which the heating of the plasma in the broken 
phase is typically lower [22]. It is possible to perform the analysis of the viability of the 
supersonic EW baryogenesis mechanism in the case of hybrid solutions, but this is left for 
future work. 

As seen from Figure [TOl in the present scenario both the f ^-ranges for which T_ > and 
for which T > become larger for increasing Higgs mass rrih- Since increasing Higgs mass 
weakens the phase transition, this might seem counterintuitive. However, for a fixed rrih, the 
^'^^-ranges become larger as ( increases, also meaning that the phase transition gets stronger. 
We conclude that the values of ruh and (, being independent parameters, have a strong 
influence on both the strength of the phase transition (as measured by Rc = h^{Tc)/Tc) and 
on the f^-range for which T > 0; the interplay among these two parameter is complicated 
and it is possible to find values of rrih and ( with the same Rc but very different (T > 0, 
f^)-range, as shown in Figure [TT] (down-left plot). 

This effect is similar to the one encountered in the analysis of section [3] for the case of 
the Bag EoS, where we saw that the (T_ > Tc)-region becomes larger as a_/a+ decrases (for 
stronger phase transitions), but this region also becomes smaller as increases (meaning 
also for stronger phase transitions). In the present case, ajv and a_ both depend on ( and 
are not independent (see Figure [TTj) . resulting in the (T_ > Tc)-region becoming larger as ( 
increases and the phase transition gets stronger. However, rrih and ( are indeed independent, 
and the relation between the strength of the phase transition and the size of the (T > 0, 
fu;)-range happens to be opposite for them. In Figure [TT] we show the dependence of various 
parameters of the phase transition (Tc, T^, T*, a^, Rc and a^/a+) with the coupling ( for 
different choices of rrih. 

Figure [12] shows the absolute value of the maximu n0 possible relative velocity between 
the plasma and the symmetric bubble wall, corresponding to ^ = Cg, in the (C, 'y«,)-plane for 

Hybrid solutions are usually allowed together with detonations [35], even though typically the pure 
detonation solution is selected by the system as its stationary expansion mode |?ni31| . However, since in the 
region below the dashed- red line in Figure [10] a detonation is not stable, the hybrid solution may be realized 
in this case. 

^^The average relative velocity is typically a factor 2 — 3 smaller, see Figure [9] 
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the particular case of nih = 200 GeV. We see that typically |(5f(cs)| ~ 10"'^ — 10~^ ^ c^, 
and diffusion is indeed efficient in the present scenario. As already discussed in section [6l 
we expect this to be a general feature for any model leading to supersonic bubble expansion 
and inverse tunneling. 
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Figure 11: Dependence of various parameters of the phase transition with the coupling ^ for 
different values of m^. Up-left: T/v (dashed line), (solid line) and (dotted line). Up-right: 
Oat. Down-left: Rc = hj^[Tc) /Tc- Down-right: a_/a+. For each value of rrih, the range of C, for 
which T, Oat, Rc and a_/a+ are plotted corresponds to the one chosen for the (^-axis in Figure [TOl 

Before concluding this section, let us briefly comment on the fact that part of the range 
for rrih analyzed here (m^ e [140 GeV - 200 GeV]) has been recently excluded by LHC for 
the case of a Standard Model Higgs [i5lH6] . In singlet scalar field extensions of the Standard 
Model it is possible to evade this constraint if the singlet scalars are light enough to open 
the decay channel h ^ Si Si. In the present case, due to the fact that the couplings of the 
scalar singlets were chosen to be universal d = (, the singlets are not light enough, since 
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^1 = C"^^ s-iid ( has to be relatively large for the phase transition to be sufficiently strong. 
However, the coupling universality = ^ is just a simplifying assumption, and it may be 
dropped (keeping ^^iVjO = ^sC) without changing the results of this section. Then, the 
constraint is easily evaded if at least one of the couplings d satisfies 



0< 



ml 



(59) 



Finally, in the general case where both the Higgs field and the singlet (s) change their 
VEVs during the transition (and after including all the relevant contributions to the free 
energy J^, such as the resummation of daisy diagrams), it is found that the phase transition 
is generically stronger [39|B0] than in the present scenario. We therefore expect that the 
supersonic EW baryogenesis mechanism may also work in the more general case, although 
a more detailed study is needed. 



Hydrodynamic Obstruction 
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Figure 12: The relative velocity between the symmetric bubble wall and the surrounding plasma 
6v as a function of the bubble wall velocity Vw for different values of (, and nih = 200 GeV. The 
solid-blue and dashed-red lines correspond to the same lines in Figure [101 namely T = and 
T_ = (they therefore delimitate the region in which baryogenesis can occur). 



8 Conclusions. 

The expansion velocity of the broken phase bubbles during a ffist order EW phase transition 
plays a fundamental role both for EW baryogenesis and for the production of a stochastic 
background of gravitational waves. Whereas fast bubble walls are needed to produce a grav- 
itational wave background detectable by future space-based interferometers, the standard 
EW baryogenesis mechanism is efficient only for scenarios in which the wall velocity is sub- 
sonic and small, since otherwise the diffusion of particle asymmetries in front of the wall is 
strongly suppressed. So far, an explicit computation of the wall velocity during a ffist order 
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EW phase transition has been performed only in the case of the Standard Model (leading 
to ~ 0.35 — 0.45 and for very small Higgs masses, already excluded by LEP), and in the 
case of a few specific models corresponding to weakly first order phase transitions (leading as 
well to rather low values for the wall velocity, Vu, ~ 0.05 — 0.1 for the MSSM). Nevertheless, 
for scenarios where the EW phase transition would be more strongly first order, qualitative 
arguments already show that the wall velocity may be much larger, thus potentially render- 
ing the standard EW baryogenesis mechanism inefficient in those scenarios, even if the three 
Sakharov conditions for baryogenesis are satisfied. 

Here we have presented a new mechanism for baryogenesis which operates for supersonic 
bubble walls, making it possible to generate a sizable baryon asymmetry when the wall 
velocity of the expanding bubbles is large. The mechanism is based on the fact that, when 
broken phase bubbles expand as supersonic detonations, the plasma in the broken phase gets 
heated up in a small region behind the bubble wall. If the temperature in this region exceeds 
the critical temperature Tc, bubbles of the symmetric phase can nucleate and grow inside 
the broken phase bubble, sphalcrons being reactivated inside these symmetric bubbles due 
to the restoration of EW symmetry. Since the symmetric bubbles are nucleated at rest in 
the moving broken phase plasma, they move through the detonation wave as they expand, 
entering at some point the region where the heating is not very efficient and the temperature 
is lower than Tg. In this region, the expansion of the symmetric bubbles is not any longer 
energetically favourable, and the bubbles start to shrink, the plasma then fiowing from 
the symmetric to the broken phase across the symmetric bubble wall. Baryogenesis then 
occurs inside the symmetric bubbles via the usual diffusion mechanism (from standard EW 
baryogenesis), and the generated baryon asymmetry is naturally transferred to the broken 
phase as the symmetric bubble shrinks. We have shown that this process happens at a 
sufficiently low velocity to guarantee that diffusion is indeed effective. 

For this supersonic EW baryogenesis mechanism to work, one further requirement is 
that a large enough volume of the plasma goes through the symmetric phase. In order to 
see whether this can indeed be the case, we have chosen a specific model where the EW 
phase transition is strongly first order (the Standard Model coupled to several singlet scalar 
fields via the Higgs portal). We have obtained the region in parameter space (in terms of the 
universal coupling ( between the Higgs and the singlet fields, the wall velocity Vyj and the 
Higgs mass rrih) where the filling factor T is large enough, and EW baryogenesis can occur 
for supersonic broken phase bubble walls. We have found that the region in which the filling 
factor is large (T — > 1) becomes smaller as nih decreases, while for a fixed it becomes 
larger for stronger phase transitions. 

The process presented here can therefore lead to significant baryogenesis when the wall 
velocity is supersonic, which may be the case in models where the EW phase transition 
is strongly first order. It could also accommodate a sizable GW signal with effective EW 
baryogenesis, leading in principle to a scenario which can explain the baryon asymmetry of 
the universe and at the same time be tested with GW detectors. 
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A. Matching Equations for Symmetric Bubble Walls 

In this Appendix we briefly analyse the matching equations across the bubble wall (jl]) for 
a symmetric bubble in the general case where the pressure is given hj p = —J^{(j),T), and 
describe the relation between the latent heat and the change in the relativistic degrees of 
freedom. We start by decomposing J-" into 

^(0,r) = ^(o,r) + J(0,T), (A.i) 

where J-" automatically fulfills ^-"(0, T) = 0. Furthermore, we can define J^(0, T = 0) = e and 
/(T) = J^(0, T) - J'(0, T = 0) = -(a/3) T^ so that 

^(0,T) = 6-^r^+J(0,T), (A.2) 

where a is related to the number of relativistic degrees of freedom of the plasma in the 
symmetric phase (it corresponds to the previously defined a+). Then, the pressures and 
energy densities in the symmetric and broken phases {psiPbi^s and Cb) can be written as 



es = aTt + e (A.3) 



p, = -T^-e-:F{ct^{T,),T,) 



dT 

We now define 



n 



A(T) 



a{T) 
1 



J(0,T)-T 



-'-b 



dT 



rp4 



(A.4) 



ce{T) is the ratio between the vacuum energy difference (between the broken and symmetric 
minima) and the thermal energy of the plasma in the symmetric phase, and A(T) is the 
liberated latent heat. For the case of symmetric bubbles the matching equations across the 
bubble wall (HI) can be written as 



VbVs 



1 - r (1 - 3a{Tb)) Vb 
3-3r(l + A(T5)) Vs 



3 + r (l-3a(Tb)) 
l + 3r (l + A(Tfe)) 



(A.5) 
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For the nucleation of symmetric bubbles, T > T^. and so a(T) as defined in (lA.4p is positive 



F (0(T), T) - (0, T) > ^ a{T) > (A.6) 

However, since the variation of J-" is the same as for the usual symmetric-to-broken phase 
tunneling (the fact that now tunneling happens in the opposite direction does not change 
the variation of J-" with temperature), we actually have 

d[T{<f>{T),T)-T{o,T)] J(0(r),r)-J(o,T) 

dT T ^ ■ ^ 

which means that, as opposed to the usual case of broken bubbles, for symmetric bubbles the 
latent heat A(T) is negative. This is however related to the fact that, when converting broken 
to symmetric phase, the number of relativistic degrees of freedom increases (increasing the 
thermal energy of the system). Comparing ( ]A.3ll with the equivalent expressions for the Bag 
EoS (dll]) and (^), we get 



eb = aT,^ + e + aT^ X{n) = a i 1 + —j + X{n) = a. (A.8) 



and it follows that in order to have a_ < a = a+, it is needed that A(Tfe) < 0. At the 
same time, the increase in the number of relativistic degrees of freedom favors the bubble 
expansion, since it contributes to the net pressure on the wall. 



B. Symmetric Bubble Growth and Stq 

As explained in section |5l due to the detonation wave's inhomogeneous background the 
symmetric bubble will be deformed from its initial spherical configuration as it expands. 
Neglecting the curvature of the broken bubble, the evolution of the symmetric bubble can 
be described by the variation of tq and zq (see Figure [13]) 



dre_ v{^)+5v{^)S^{e) dze_ (ve^^ 



l + v{'f) 5v (^) Sin(^) 

Sin(g) = ~ ^'^ for e>0 Sin(g) = - ~ for ^<0. (B.2) 

2 — 

where Sv{C,) is given by (j33D- Also, from Figure US] one has = and rn = rout as 

2 2 

defined in ( |35|) . and it is then clear that (!35|) derives from (]B.1|) . 



Let us now turn to the evaluation of the volume of plasma that goes through a symmet- 
ric bubble. Focusing on the broken bubble's radial direction, this amounts to obtain the 
quantities dr^^^ and Sr^^^ as a function of tq and ^ (as defined in fl37p ). This can be done 
using the self-similarity of the detonation wave. We start by defining as (see Figure [T^ 

rm(ro,Oi) = ^c(ro,Oi) • (B.3) 
Then, due to the self-similarity of the detonation wave, and using ( !T0|) and ( ]B.3p . the ratio 
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Figure 13: Parametrization of the symmetric bubble evolution. 



does not depend on the value of tq (see Figure [14]). Since C(^/J — 1 is already independent 
of To, we get 



Using similar arguments, one also has 



To 



(ro— independent) . 



(ro— independent) . 



Then, from f ITOl) and f l371) we obtain 



ro 1 



Sr. 



\ C({) 



Finally, we arrive at 



5ro = Srl^^ + 6rl^^ 



c{i) 



(B.5) 



(B.6) 



(B.7) 



(B.8) 



As a final remark, even though fIB.ip can be used to study the deviation from spherical 
symmetry of the symmetric bubbles as they expand, for all practical purposes concerning 
the mechanism of EW baryogenesis outlined in section [H] one can neglect its effects on the 
computation of the volume of plasma that goes through a symmetric bubble Vs and assume 



Vs ~ {Srof 



(B.9) 
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t 



Figure 14: Detail of Figure [U where the scale-invariant (self-similar) ratios used to prove ()10p and 
the ro-independence of (|B.5P and ()B.6P are shown explicitly. 
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